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Abstract 

The equation of state for neutron matter is calculated within relativistic Hartree-Fock approxi- 
mation. The tensor couphngs of vector mesons to baryons are included, and the change of baryon 
internal structure in matter is also considered using the quark-meson coupling model. We obtain 
the maximum neutron-star mass of ~ 2.0Mq, which is consistent with the recently observed, pre- 
cise mass, 1.97ibO.O4M0. The Fock contribution is very important and, in particular, the inclusion 
of tensor coupling is vital to obtain such large mass. The baryon structure variation in matter also 
enhances the mass of a neutron star. 

Keywords: neutron stars, equation of state, quark-meson coupling model, relativistic Hartree-Fock calcula- 
tion, tensor interaction 



Electronic address: koichi.saito@rs.tus.ac.jp 



Because the internal structure and the mass of a neutron star significantly depend on the 
equation of state (EOS) for neutron matter, the observation of pulsars provides important 
information to understand the properties of dense matter 11]. At sufficiently high density, 
nuclear or neutron matter would be composed of not only nucleons (A^) and leptons {£) 
but also hyperons (Y) jl-5| and, possibly, meson condensates {2, 1^- Since it is important 
to calculate the EOS relativistically, many studies on neutron matter have been performed 
within relativistic Hartree approximation. 

Recently, the precise observation finds the pulsar with the mass of 1.97 ± 0.04Mq (solar 
mass) It is, however, difficult to explain the observed mass by the EOS calculated in 
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and thus reduces the maximum mass of a neutron star [3|-|5|. Thus, it seems very urgent 
to consider how this discrepancy can be reconciled. It may be necessary to first study the 
effects of the Fock term, the tensor couplings of vector mesons and the baryon structure 
change in matter, and see how those effects contribute to the EOS or the maximum mass of 
a neutron star. 

In fact, the relativistic Hartree-Fock calculation of the ^OS for neutron matter has already 
been performed in Quantum Hadrodynamics (QHD) [8] or in the quark-meson coupling 
(QMC) model 9|, and it is pointed out that the Fock term is significant to enhance the 
maximum mass. However, in Ref . |8| , the effect of the baryon structure change in matter is 
ignored. In Ref. jol, although the maximum mass of about 2Mq is obtained, their calculation 
misses the tensor coupling and the space component of baryon self-energy. 

In this paper, we study the properties of nuclear and neutron matter within relativistic 



Hartree-Fock approximation 10|, lll| , where the pion exchange and the tensor couplings of 



vector mesons to baryons are included in the Fock term. Furthermore, using the QMC 
model 12| , we consider the variation of the quark substructure of baryon in matter. In the 



past few decades, the QMC model has been extensively developed and applied to various 
nuclear phenomena with tremendous success 12|]. We also apply the chiral QMC (CQMC) 
lich is recently extended to include the quark-quark hyperfine interaction 



model [13|,Q, w 
due to the gluon 
matter. 



15| and pion exchanges based on chiral symmetry, to the EOS for neutron 
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The Lagrangian density for dense matter is given by 



(1) 



where 

I^B = Y, ^B{il^.d^ - MB)i^B, C,i = Y. M^f.d'^ - (2) 

B I 

with ipBi^) the baryon (lepton) field and MB{mf) the free baryon (lepton) mass. The sum B 
runs over the octet baryons, p, n, A, S^*^^ and S*^", and the sum i is for the leptons, e~ and 
fi-. For the baryon masses, we take = 939 MeV, Ma = 1116 MeV, Ms = 1193 MeV 
and M= = 1313 MeV, respectively. 
The meson term reads 

Cm = \ - mla^) + ^m^^o;^ - \w^,W^- 

+ Imjp, ■ r - \r,u -R^' + l (9,7f ■ a^Tf - mln') , (3) 

with 

Wf,u = d^Uu - dyUf,, R^u = d^pu - d^Pf,, (4) 

where we consider the isoscalar {a and u) mesons and the isovector (tt and p) mesons, and 
the meson masses are respectively chosen as = 550 MeV, m^^ = 783 MeV, = 138 
MeV and nip = 770 MeV. 



The interaction Lagrangian is given by 



^B 

B 



2M 

^ ^B, (5) 



where Ib is the isospin matrix for baryon B (we set Ib = when B is iso-singlet) and the 



common, scale mass Ai is taken to be the free proton mass [16|. The cr-, oo-, p-, ir-B coupling 
constants are respectively denoted by g^jBi^r), QujB, QpB and /ttb, while f^iB and fpB are the 
isoscalar- and isovector-tensor coupling constants. 

In QHD, QfjB is constant at any density because the baryons are assumed to be structure- 
less. In contrast, in the QMC and CQMC models, the ganio') has the a-field dependence. 



which reflects the baryon structure variation in matter 



In relativistic Hartree-Fock 



approximation, the a-field is replaced by the constant mean-field value, a. Furthermore, for 



TABLE I: Values of and for various baryons in the QMC or CQMC model. 





QMC 
as (fm) 


bB 


CQMC 

as (fm) 


bB 


N 


0.179 


1.00 


0.118 


1.04 


A 


0.172 


1.00 


0.122 


1.09 




0.177 


1.00 


0.184 


1.02 




0.166 


1.00 


0.181 


1.15 



simplicity, we here use the parameterization of the quark scalar-density ratio, Cb{o'), in the 
linear form [l2, 



Cij(o-) = 1 - Ob X (g^NO-), 



(6) 



where g^N is the a-N coupling constant at zero density and is a parameter. Using this 
parameterization, the a-B coupling constant in Eq. ([5]) is given by 12, llTI] 



gaB{(^) = gaBbs 1 - ^{gaN(^) 



(7) 



where we introduce a parameter 6^ to use Eq. ([7]) in both the QMC and CQMC models. 
The values of and bs are listed in Table [H If we set = and 6^ = 1, gasip^) is identical 
to the a-B coupling constant in QHD. 

In the present calculation, we add a non-linear (NL) term of the a. 



TTi-\ 92-3 . 5^3-4 



(8) 



to only the QHD Lagrangian, because the EOS given by the naive QHD is too hard 10 |. 
We call this QHD+NL. 

To sum up all orders of the tadpole and exchange diagrams in the baryon Green's function, 
Gfi, we use the Dyson's equation 



GB{k) = Gl{k) + Gl{k)T.B{k)GB{k), 



(9) 



where is the baryon self-energy and G% is the Green's function for the free baryon. In 
nuclear or neutron matter, the baryon self-energy is generally written as 



T.B{k) = J:%{k) - 7oS°j(A;) + (7 ■ k)El{k), 



(10) 



TABLE II: Functions Ai, Bi, Ci, and Dj. The index i is specified in the left column, where 
V{T) stands for the vector (tensor) couphng at each meson-BB' vertex. The last row is for the 
(pseudovector) pion contribution. 



Di 



a 

UJVV 
UJnpT 
OJVT 
PVV 

Ptt 

PVT 



2 „ 



fuiB 

2M 



2M j '"'P^P 



( 2M ) "^a)®<^ 



-Hb^p 
-3 llf^Ymle 



2M 



p^P 



A ( ikg. ) 11/ 



-^g^B^P 



a(IpB. 
^\2M 



fnB 



1 O fuiBdujB p 
2M 



1 r, fpBgpB 

2M ^ P 



with k the unit vector along the momentum k and S^^°^'^' the scalar part (time component) 
[space eompouenti of the self-euergy. Furthermore the effective baryon ,„as. „„,„e,.um 

and energy in matter are respectively defined by [UJ 



Mm 



h.*p 

^B 



MB + ^%{k), 

(U*0 t* \ _ ( 
y^B ! '^B) 



(A;" + E°5(fc),A; + A;S^(fc)), 

1/2 



The baryon self-energies in Eq. ( JTOl) are then calculated by 11| 



E%{k) 



-gaB[(T)Cr 



'■BB' 



B' 



7l Jo 



kp , 

" R' 



dqq 



(11) 
(12) 

(13) 



MMsls-fk a) + ^—D ia k) 



Ss(^) = -9ujBUJ - gpB{lB)3p 



B'i 



BB' 

7^ (4vr)2fc Jo 



dqqAi{k,q), 



7l (4vr)2A; ^ 



(14) 
(15) 
(16) 



where Co and p are respectively the mean-field values of u and p mesons in matter, kpg is 
the Fermi momentum for baryon B and the factor, {Ibb'Y^^^ isospin weight at the 

meson-i?i?' vertex in the Fock diagram. The index i in the sum and the functions Ai, B^, 
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Ci and Di in Eqs. (fT^ - (fT6|) are explicitly given in Table HTj in which the following functions 
are used 

"mf + {k + qY' 



&i{k,q) = In 



Akq 



mf + {k- 
(e + q^ + m^i) ei{k,q)-l, 



%{k,q) = [{k^ + q^-m^/2)^,{k,q)-kqe,{k,q)], 
n,(fc,g) = [{k'' + q^)^,{k,q)-kqQ,{k,q)], 
Ti{k,q) = [kQ,{k,q)~2q^,{k,q)]. 

The mean-field values of u and p are respectively given by the usual forms 



(17) 

(18) 
(19) 
(20) 
(21) 



— Pb, P=V— (^b)3Pb, 22) 
m^, ^ — ^ m, , 

where the density of baryon B is ps = (^^ef^-'^^s ^i^h the spin of B. 

On the other hand, combining with Eqs.(lT4l)- (fT6ll . the value of a is self-consistently 
calculated by [8] 

^ = E ':SbBCBia)p% - ^{g,-a^ + g,a% (23) 



where the scalar density of baryon B reads 

, 2Jb + 1 



Pb 



IB 

27r2 



,M^(/.) 



El{k) 



(24) 



In Eq.([23D, §2 and 5(3 vanish in the QMC or CQMC model, while, in QHD+NL, hsCBi^) is 
taken to be unity. 

The total energy density (pressure) for matter can be divided into the baryon and lepton 
contributions, namely e = + Q {P = Pb + Pi), where the baryon energy density, e^, and 
pressure, Pb, are expressed by 



_^ 2Jb + 1 

(2vr)3 



B 



dk 



TB{k) + U^B{k) 



o-^ (92 , 93- , 



(25) 



with 



MBMUk) + kki MUk)^%{k) + k^^EUk) ^ 
^^^^^ EUk) ' ^^^^^ " EUk) 



and 



2 9 { tB 



Ps = n%-—\ — 
duB KriB 



(26) 



(27) 



TABLE III: Coupling constants and calculated properties of symmetric nuclear matter at n^. In 
the left column, H (HF) stands for the Hartree (Hartree-Fock) calculation. The incompressibility, 
Ky, and the symmetry energy, 04, are in MeV. The coupling constants, (72 (in fm~^) and 5(3, in 
QHD+NL are chosen so as to adjust and to the values given by the QMC model. 



Model 








92 


93 




M*^/Mn 


04 


QHD+NL(H) 


5.34 


5.41 


1.88 


11.07 


95.59 


280 


0.80 


36.9 


QMC(H) 


5.45 


5.41 


1.88 






280 


0.80 


36.9 


CQMC(H) 


5.75 


7.11 


1.82 






302 


0.76 


36.9 


QHD+NL(HF) 


3.63 


6.66 


0.48 


16.27 


48.52 


280 


0.73 


38.6 


QMC(HF) 


3.50 


6.41 


0.48 






280 


0.73 


36.9 


CQMC(HF) 


3.59 


7.34 


0.48 






300 


0.70 


37.4 



In Eg. (1271). the total baryon density, ub, is given by ub = J2bPb- 

The numerical result for the properties of symmetric nuclear matter is presented in Ta- 
ble lllli In the present calculation, the a-N and u-N coupling constants are determined so 
as to reproduce the saturation energy (—15.7 MeV) at normal nuclear density (n^ = 0.15 
fm"^). 

In the Hartree calculation, the p-N coupling constant, gpN, is determined so as to fit the 
symmetry energy, 04 = 36.9 MeV Q, O, [isl, and the meson-hyperon coupling constants, g^Y, 
g^Y and gpY, are taken to be the values derived from the quark model [3|, 

In the Hartree-Fock calculation, we use the more reliable set of the coupling con- 
stants (except ga-B and q,.,/v) determined in the Nijmegen extended-soft-core (ESC) model 
(see Table VII in Ref. 16|]). Furthermore, from the recent analyses of hypernuclei 
and hyperon production reactions, it is inferred that the A (S) [S] feels the potential, 
f^A(s)[s] — —30 (+30) [—15] MeV, in nuclear matter 

coupling constants, g^^, g^j^ and (^o-h, so as to reproduce the suggested potential values 
at n%. In the case of QHD+NL (QMC) [CQMC], we take g„K/y/A^ = 2.88(2.78) [2.44], 
g^s/V^^ = 2.05 (1.97) [1.94] and g^s/V^ = 2.06 (1.98) [1.70]. 

In Table IIVI the contents of the nucleon self-energies in the QMC model are presented. 
The Fock term contributes to the self-energy significantly. In symmetric nuclear matter, the 
effect of the tensor coupling of u meson is small, while that of the p meson is considerable; 
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TABLE IV: Nucleon self-energies, at in symmetric nuclear matter. The values of the 

self-energy are in MeV. 





ys 




yv 


ndi Li 


- lOU 




n 

u 












-131 


-151 





Exchange o" 


13 


-14 







-63 


-32 


-3 


TT 


-4 


4 


-4 




-69 


13 


14 


total 


-255 


-180 


6 




for example, the p-meson contributions of the vector iVV), tensor (TT) and vector-tensor 
iVT) mixing to the scalar part, S^, are respectively —14 MeV, —58 MeV and 4 MeV. The 
space component, S}!^, at n% is relatively small. 

In Fig. [1], we show the nucleon self-energies calculated by the QMC model as functions 
of the total baryon density, ub- The scalar and time components, S^", in the Hartree-Fock 
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FIG. 2: Particle fractions in neutron matter. 



TABLE V: Neutron-star radius, R 

max (in km), the central density, (in fm "^j, and the ratio of 
the maximum neutron-star mass to the solar mass, Mmax/^Q- The Hartree and the Hartree-Fock 
calculation with (without) hyperons are denoted by npY (np). 







np 






npY 






Rmax 


ric 


Mmax/MQ 


Rmax 




Mmax/MQ 


QHD-FNL(H) 


11.3 


1.04 


2.00 


12.5 


0.86 


1.56 


QMC(H) 


11.5 


1.01 


2.05 


12.5 


0.86 


1.60 


CQMC(H) 


11.8 


0.92 


2.20 


12.5 


0.88 


1.66 


QHD-FNL(HF) 


11.7 


0.95 


2.15 


11.9 


0.95 


1.92 


QMC(HF) 


11.5 


0.97 


2.11 


12.0 


0.92 


1.95 


CQMC(HF) 


11.9 


0.90 


2.23 


12.3 


0.87 


2.02 



calculation are much deeper than those in the Hartree calculation. It is noticeable that the 
space component also becomes deep at high densities and thus it is not negligible in dense 
matter. 

In a neutron star, the charge neutrality and the /3 equilibrium under weak processes 
are realized. Under these conditions, we calculate the EOS for neutron matter and solve 
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FIG. 3: Neutron-star mass as a function of the radius. The left (right) panel is for the Hartree 



(Hartree-Fock) calculation. In solving the TOV equation, we use the EOS of BBP 
at very low densities. 



2m and BPS 



the Tolman-Oppenheimer-Volkoff (TOV) equation. In Fig. |2], we show the QMC result of 
particle fractions in relativistic Hartree-Fock approximation. With respect to hyperons, only 
the S~ appears around 0.49 fm~^ and the other hyperons are not produced at densities below 
1.2 fm~'^. We note that, in the case of QHD-I-NL, the S~ first appears around 0.43 fm~^, 
and that the A and H° are produced at densities beyond 0.69 fm~'^. This tendency is similar 
to the result of Ref. It may be interesting to compare the present result with the Hartree 
result given in Ref. 19 |. 

We summarize the properties of neutron star in Table |Vl and show the mass as a function 
of the neutron-star radius in Fig. [3l As known well, the inclusion of hyperons generally 
reduces the mass of a neutron star. However, because the Fock contribution makes the EOS 
hard, the maximum mass in the present calculation can reach the recently observed value, 
1.97±0.04Mq. If we ignore the tensor coupling in the Fock term, the difference between the 
maximum masses in the Hartree and the Hartree-Fock calculations is not large. Therefore, 
the tensor coupling (especially, in the high density region) is very vital to obtain the large 
neutron-star mass. The variation of the quark substructure of baryon in matter also enhances 
the mass, because the QMC model includes the (repulsive) many-body interaction given 
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through the scalar polarizabihty in the couphng constant, gaB{(^) 0, Q|- Furthermore, the 
quark-quark hyperfine interaction due to the gluon and pion exchanges increases the mass 
of a neutron star. We should note that the inclusion of the meson-baryon form factor at the 
interaction vertex in the Fock term does not change the present result much. 

In summary, we have studied the effects of the Fock term, the tensor couplings of vector 
mesons and the baryon structure variation on the properties of a neutron star. The present 
calculation can produce the maximum neutron-star mass of ~ 2.0Mq, which is consistent 
with the recently observed mass, 1.97 ± O.O4M0. The Fock contribution is very important 
and, particularly, the inclusion of tensor coupling is necessary to obtain such large mass. 
The in-medium variation of baryon structure also makes the EOS hard and thus it enhances 
the mass of a neutron star. In the future, it may be desirable to consider the degrees 
of freedom of A isobar [l|, O], the i^'- meson exchange and the baryon mixing in the Fock 
diagram. It is also interesting to study the possibility of meson condensates [l|, l2|, |6[ and/or 
quark matter 

The authors would like to thank K. Nakazato, C. Ishizuka and H. Suzuki for fruitful 
discussions on the equation of state for a neutron star. 
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